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The structure of a weak shock in a quantum plasma is studied, taking into account 
both dissipation terms due to thermal conduction and dispersive quantum terms due 
to the Bohm potential. Unlike quantum systems without dissipations, even a small 
thermal conduction may lead to a stationary shock structure. In the limit of zero 
quantum effects, the monotonic Burgers solution for the weak shock is recovered. 
Still, even small quantum terms make the structure non-monotonic with the shock 
driving a train of oscillations into the initial plasma. The oscillations propagate 
together with the shock. The oscillations become stronger as the role of Bohm 
potential increases in comparison with thermal conduction. The results could be of 
importance for laser-plasma interactions, such as inertial confinement fusion plasmas, 
and in astrophysical environments, as well as in condensed matter systems. 

I. INTRODUCTION 



Quantum plasmas, where the finite width of the electron wave functions gives rise to 
collective effects jl, 2, , are currently a rapidly growing field of research. Many of the studies 
are motivated by the potential for application to nanoscale systems |4|], such as quantum 



wells M, u 

n 

sources |9|, 



7J , laser fusion plasmas 



next-generation high intensity light 



tracold plasmas 

lo| , and plasmonic devices . Moreover, nonlinear effects in quantum plasmas, 



such as the formation of dark solitons and vortices [12j, interaction between quantum plasma 
osculations and electromagnetic waves . quantum turbulence [ijj], and solitary structures 



TgI supported by the electron spin 



17 



181 ]. are currently in focus as well. 



There has also been much interest in shocks in quan tum-like systems, such as nonlinear 
optical fibers and Bose-Einstein condensates 



19 



20, 



21 



22] . The structure of such quantum 



shocks is markedly different from the classical ones |23[. The shock structure in classical 
fluids/gases is governed by transport processes, i.e., the viscosity and thermal conduction. 
A classical shock propagating with constant velocity displays a stationary structure. If the 
shock is weak, then transition from initial matter to compressed one may be described by 



2 



the smooth monotonic Burgers solution 23|]. In contrast to classical fluids, quantum media 



typically exhibit dispersion due to the Bohm potential instead of dissipation 



24 



25 



19 



20 



21 



22 



27 



281 ] . For this reason, even a quantum shock propagating with constant velocity 



in a uniform medium does not posses a stationary structure. Transition from initial to 
compressed quantum media happens in the form of a train of solitons of different amplitudes 



19 



20 



21 



22J. The solitons propagate with different velocities, which makes the whole 



structure intrinsically non- stationary. Obviously, a train of solitons also provides a non- 
monotonic transition from initial to final state of the medium. However, there are quantum 
systems with both dissipations and dispersion, such as quantum plasmas. The viscosity in 
plasma is determined by ions and it is typically negligible. Still, electron thermal conduction 
may be quite strong both in classical and quantum plasmas 25[]. Therefore, shocks in such 
plasmas may demonstrate transitional behavior between the classical and quantum domains. 
The purpose of the present paper is to trace such a transition by studying weak shocks. 

Here, we derive a nonlinear equation governing the structure of a weak shock in quantum 
plasma. The equation contains both dissipation terms (due to thermal conduction) and 
dispersive quantum terms (due to the Bohm potential). Unlike quantum systems without 
dissipation, even relatively weak thermal conduction may lead to a stationary structure of a 
shock. In the limit of zero quantum effects we recover the monotonic Burgers solution for the 
shock structure. Still, even small quantum terms make the transition non-monotonic with 
the shock driving a train of oscillations into the initial plasma. The oscillations propagate 
together with the shock with the same velocity. The oscillations become stronger as the 
role of Bohm potential increases in comparison with thermal conduction. The oscillations 
resemble the soliton train in quantum shocks without dissipations. 



II. GOVERNING EQUATIONS 

The basic equation of nonrelativistic quantum mechanics is the Schrodinger equation. The 
dynamics of an electron, represented by its wave function ip, in an external electromagnetic 
field (cf), A) is governed by 

^4 + 4~ (V + %A)\ + eU = 0, (1) 



dt 2m P \ h 
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where % is Planck's constant, m e is the electron mass, and e is the magnitude of the 
electron charge. This complex equation may be written as two real equations, writing 
ip = y/p expiS/h, where p is the amplitude and S the phase of the wave function, re- 
spectively 26(| . Such a decomposition was presented by de Broglie and Bohm in order to 
understand the dynamics of the electron wave packet in terms of classical variables. In Ref. 



251 ] the Wigner function was employed for the purpose of obtaining a set of quantum hydro- 



dynamic equations. In this way, an arbitrary number of conservation equations, in particular 
an energy conservation equation, may be obtained before closure. Here we will just briefly 
review the Bohm-de Broglie approach, making use of the energy conservation equation from 
Ref. (25( 1 . Using the decomposition of the wave function in terms of its amplitude and phase, 
Eq. (JH) gives 

| + V-(pu)=0, (2) 

and 

me _ = e( E + uxB) + — V(^), (3) 

where the velocity is defined by u = VS/m e , and E = — V0 — dtA and B = V x A. The 
last term of Eq. ([3]) is the gradient of the Bohm-de Broglie potential, and is due to the effect 
of wave function dispersion. We also note the striking resemblance of Eqs. (j2j) and ([3]) to 
the classical fluid equations. 

Suppose that we have N electron wavefunctions, independent apart from their interaction 
via the electromagnetic field. For each wave function ifj a , we have a corresponding probability 
V a . From this, we first define ip a = y/~p^exp(iS a /h) and follow the steps leading to Eqs. (J2J) 
and (131). We now have N such equations the wave functions {i/j a }- Defining 

TV 

p=J2 V aPa (4) 

a=l 

and 

N V o u 

U = (Ua) = 2^ " , (5) 

a=l P 

we can define the deviation from the mean flow according to 

w Q = u a - u. (6) 

Taking the average, as defined by (jSJ), of Eqs. and and using the above variables, we 
obtain the quantum fluid equation 

^ + V-(pu) = (7) 



4 



and 



Pe 



' d_ 

dt 



u R ■ V u R 



ep e 



fE + u e x B) - VP e + 



h 2 p c 



V 



'Po 



2m e \ V \JPol ) I 

where we have assumed that the average produces an isotropic pressure P = p e (|w Q | 2 ) We 
note that the above equations still contain an explicit sum over the electron wave functions. 



For typical scale lengths larger than the 
term by the Bohm-de Broglie potential 

/V 2 , 



ermi wavelength Xp, we may approximate the last 



25] 



'Po 



V 2 Jp e 



(9) 



\ y/P^ I "i\fP~e 

where the factor 1/3 comes from an isotropic averaging. 

Treating the ions as fully classical due to their large inertia, we can derive a a set of single- 
fluid quantum magnetohydrodynamic (MHD) equations, following [18] where a similar set 
of equations for a perfect fluid type quantum MHD system was derived. Thus, we obtain 



the equations of mass, momentum and energy transfer in a quantum plasma 18 

0. 



2a 



dp d 



-(pu,) + —(pu^) 



dP 



+ 



ti 2 



a 



d 2 



dxj 12m e rrii dxi \ dxjdx, 



In p 



(10) 



(11) 



and 



d_ 
dt 



u 

Pe + P- 



-pui 



ti 2 



pV 2 In p 



24m e m L 
h 2 d 2 \np 
12m e mj dxjdxi dx 



d 




dxj 


pUj ^ 



u 

T 



ti 2 



2Am f ,m i 



V 2 lnp 



dT 



j J 



0. 



(12) 



where p is the fluid mass density, Uj is the fluid velocity, P is the pressure, e, h are thermal 
energy and enthalpy, k is thermal conduction, and m.j, m e are the ion and electron masses. 
The energy conservation equation ([12"]) is derived using the Wigner approach 25| for the 
electron dynamics and combining this with the ion equation in the MHD limit. Here we 
have neglected the effects due to the magnetic field, that are assumed small in comparison 



The hydrodynamic 



to the other governing terms. Such terms can easily be included [181 , 
equations should be complemented by the thermodynamic equation of state. We take the 
equation of state to be that of an ideal gas, i.e., 

r, 7-1 



7 



-CppT, 



(13) 



5 



and 

h = C P T, (14) 

and let the electron thermal conduction k be oc T 5 / 2 . Here Cp is heat capacity at constant 
pressure and 7 is the adiabatic exponent. We stress that the forms ( 1T31) and (1T41) play a 
minor role for weak shocks. 



III. SHOCK SOLUTIONS 

We consider a planar stationary shock. In the reference frame of the shock, Eqs. (fT0l) - (|T2l) 
may be integrated as 

pu = p u , (15) 
P + -12^"V = P » + ^ (16) 

and 

m 2 ft 2 d 2 lnp k dT _ h u o (17) 
2 8m e mi dx 2 pou dx 2 ' 

where the subscript refers to the uniform plasma ahead of the shock and uq is the shock 
speed. As we can see from ([I5]) - ([PT]) . quantum effects do not influence the properties of the 
uniform flow behind the shock; they are important only for the shock structure. Next, we 
introduce the parameters 

CpPqUq ^ 

and 

Ma 2 = ^0, (19) 

and the scaled variables p/p = u /u = R, T/T = 1 + r\ = x/L. Here L is the 
characteristic length scale determined by thermal conduction; in the classical case it may be 
treated as the shock width with the accuracy of a numerical factor of order unity. The other 
parameter is the Mach number, Ma, which compares the shock velocity u$ to the initial 



sound speed y^/Po/po in the plasma and characterizes the shock strength. The parameter 
d denotes the deviation of the temperature, produced by the shock wave, from the initial 
value. In the case of weak shocks we have Ma — 1 <C 1, ie., the shock velocity marginally 
exceeds the sound speed, and as does the temperature from the initial value Using 



6 



the scaled variables, we reduce Eqs. (IToT) - (ITTj) to 



s 7Ma 2 „ d 2 lnR „ 
R(l +$) + - Q R ^T = 1 + 7Ma 2 , 



and 



tf + ^Ma 2 -(l + tf) 5 / 2 ^-^ 
2R 2 K J drj 2 7 



1) d 2 \nR 



7 



dr] 2 



-Ma 2 



where 



Q 



h 2 Po L 2 



(20) 



(21) 



(22) 



12m e rriiPo 

is the parameter comparing the role of quantum and classical effects in the shock dynamics. 
This parameter can be interpreted as a quantum Mach number. The system ( 120]) - (I2~T]) 
determines the structure of a shock wave in a quantum plasma. In the classical case we have 
Q = 0, and Eq. (T2"0~l) gives an algebraic relation between the density and temperature 



1 

R 



1 + 7Ma 2 
2 7 Ma 2 



1 ± 



\ 



1 



4 7 Ma 2 (l + tf) 



;i + 7 Ma 2 ) 



2^2 



(23) 



The positive sign in (1231) gives rise to shock solutions, while the negative sign corresponds 
to deflagrations [23(. We note that the density and temperature of the compressed matter 
increase together in a shock. In deflagrations, the temperature increase leads to a decrease 



of the density, as in, e.g., laser ablation and flames [30j, |31|, |32j, l33j. Substituting (12"3"|) into 
( f2T|) . we obtain a single differential equation for the temperature in a classical shock. In the 
case of strong quantum shocks, one has to solve a system of two differential equations. 



A. Weak shocks 

In the present paper we investigate only the case of a weak shock with Ma 2 - l = /i<l 
and t) < 1. A more general case will be studied elsewhere. As note above, this value of 
\i characterizes a shock velocity marginally above the sound speed. In the case of a weak 
shock in the linear approximation, Eq. (123]) may be simplified according to 

1 i? 

- = 1 \ . (24) 

R 7 Ma 2 -l V ; 

Taking into account the quantum dispersion and the weak nonlinearity in (f20|) . we find 

R~ ~ 7 Ma 2 - 1 + 7 Ma 2 - 1 drj 2 ( 7 Ma 2 - l) 3 ' ^ ' 
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Equation (125]) relates the density to the temperature in a weak shock. Taking into account 
the linear approximation (|24|) . we may simplify the quantum dispersive term as 

dHnR 1 d 2 d 



drf 7 Ma 2 - 1 drf 

Substituting ([25]) into (J2U), we find 



(26) 



/ 7 + lV 2 3- 7 ^d 2 tf . ^dd 

^-{^)Y + ^r Q w = h - 1) cT v (27) 

Equation ( |27j) describes the structure of a weak shock in a quantum plasma with finite 
thermal conduction. Dissipation and quantum terms enter Eq. (T2"7l) as the first and second 
derivatives of the scaled temperature. All derivatives tend to zero in the uniform flows 
corresponding to the initial and final plasma states at rj — > ±00. Taking into account that 
$ = in the initial flow, we find the relation between scaled temperature increase in the 
shock wave and the scaled shock speed 

Thus, we can rewrite (l2Tj) as 

W . fl= M* HH£ (29) 
v ' 7 + 1 dr) 7(7 + 1) dr] 2 K ' 

The solution to fl29l) changes from $ = at 77 — > —00 in the initial plasma ahead of the 
shock to •& = $1 at rj —>■ 00 in the compressed plasma behind the shock. 



B. Classical/quantum transition in the schock 

We are interested in solution to (129]) for any parameter value Q from (classical plasma) 
to infinity (quantum plasma without dissipations). To simplify our study of Eq. (129]) . we 
may rescale the temperature according to <f> — "d/fii, so that changes from in the initial 
plasma to 1 in the compressed plasma. We also rescale the coordinate 

^ _r? 2(7-l) 2 ' [ } 



and introduce a new parameter 



(7 + 1X3-7) 

g = q * 27(7-1)3 » ( 31 ) 
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Then f|29l) reduces to a concise form 

+ (32, 

The parameter g describes the relative role of quantum effects and thermal conduction in 
the shock. Equation (1321) is the main result of our paper. 

In the case of zero quantum effects (q — 0), Eq. (1321 goes over to the stationary Burgers 
equation with the solution 

6XP ^ (33) 



1 + exp £ 

The influence of quantum effects may be analyzed analytically in the limit of small q ^ 1. 
Because of the invariance to a shift in space, the solution to fl32l) may be presented the form 

= m > de = f ^ (34) 



so that 



:w)=/ (35) 



To th order in q (^C 1) we have f(6) = 0(1 — 6), and to 1 st order in q, Eq. (13"21 becomes 

6(1-6) d6 



(36) 



l-q(l-26) di 

The classical solution ( J33l) was symmetric in space with respect to the central point 6 = 1/2. 
As we can see from ( 1361) . quantum effects make the scaled temperature slope d6/d^ steeper 
in the front part of the shock, for 6 < 1/2, and smoother at the back side, for 6 > 1/2. 

We can also solve (|32|) analytically in the limit of strong quantum effects, q — > oo. In 
that case the leading terms in (132]) are 

0(1 " <P) = -9^- (37) 

Integrating ( 137|) once, we obtain 

*l^l + C= -l(f), (38) 
2 3 2 \d£) K ' 

where C = —1/6 due to the boundary condition at the back side of the shock (6 — > 1). 
Equation (|38|) may be also rewritten as 

(0-l) 2 (20+l) = 3g(^f) , (39) 



d£ 
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and integrated to give 



^tanh 2 (^)-i (40) 



The solution (1401) is a dark soliton of Korteweg-de Vries type (similar solutions have previ- 
ously been found in quantum hydrodynamics [13, [16] ). The solution fT40l) is characterized by 
a new length scale y/q; and it tends to unity — > 1, for £ — > ±oo. We note that neglecting 
dissipation, we cannot come from initial state = to the final state 0=1 smoothly 
Therefore, the shock inevitably contains a weak discontinuity, which is a surface where 
is continuous but d(j)/d^ has a discontinuity. This weak discontinuity develops at the front 
side of the shock where = and 

d ± = ± J-. (41) 
ci£ 3^q { ' 

The temperature profile may reach the point where = either from the "bright" or "dark" 
side, depending on the positive and negative sign, respectively, in Eq. (14TI) . However, in both 
these cases the transitional region is just a part of the dark soliton solution (T40|) . see Fig. [IJ 
The weak discontinuity may be removed taking into account a small but finite dissipation. 
In the limit of strong quantum effects and weak dissipation it is more convenient to rescale 
the space variable as ( = ^/^/q. In that case, the dissipation in Eq. (1321) is small (as 

i/V?<i) 



:w) = 4=£-^- (42) 



Dissipation modifies the soliton solution Eq. ( 1401) on the front side at £ —>■ — oo. Because of 
the dissipation, cannot reach unity at £ — > — oo, but tends to zero in the form of decaying 
oscillations. When is close to zero, Eq. (14"2"I) describes small linear oscillations 

1 d0 d 2 ^ 



^ Vvd( dC 

decaying at ( — > — oo according to 



(43) 



oc exp 



: + 2^) c 



(44) 



Equation (1321) is reduced to a form (1431) for any non-zero value of g as soon as the temperature 
comes sufficiently close to the initial value, 1, at the front side of the shock. Therefore, 
even for a small but non-zero quantum effects we should expect oscillations ahead of the 
shock. The oscillations decay quite fast in the case of relatively small q, but they form 
a long wave in the limit of large q. Numerical solution to Eq. (1321) is shown in Figs. [2, 
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[3] for different values of the parameter q = 0; 1; 5; 1000. The plot with q = shows the 
Burgers solution, which describes a monotonic transition from the initial to the compressed 
plasma in a classical weak shock. In the case of small but non-zero quantum effects, q = 1, 
we can see one well-pronounced "dark" region with temperature below the initial value 
(<f) < 0). Still, the oscillations decay fast for q = 1 and they may be observed only on a 
small scale of 1. Increasing the role of quantum effects, q = 5, we can clearly see a 
number of peaks and troughs ahead of the shock wave. Finally, in the case of large quantum 
effects, q = 1000, the front side of the shock looks like a train of oscillations decaying at 



( — > — oo. The 



medium 19, 



20. 



ast plo t resembles a non-stationary train of solitons in a purely quantum 



21 



221 ] . Still, we would like to stress that the shock structure obtained in 
the present paper is stationary; the train of oscillations propagate together with the shock 
with the same velocity. 

IV. SUMMARY 

In this paper, we have investigated the classical-quantum transition in weak shocks, using 
a quantum fluid model with finite heat conduction. Both analytical and numerical results 
were presented, and it was found that soliton trains can occur at the shock front in the 
quantum regime. Such significant modifications of the front structure could be of interest 
in laser fusion plasmas. 
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